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ABSTRACT 



J> ' Scattering of starlight by dust, molecules and electrons in spiral galaxies will 

fS| , produce a modification of the direct intensity and a polarization in the ob- 

. served light. We treat the case where the distribution of scatterers can be 



considered to be optically thin, and derive semi-analytic expressions for the 



' resolved intensity and polarized intensity for Thomson, Rayleigh, and more 



general scattering mechanisms. These expressions are applied to a parametric 



Q . model spiral galaxies. It is further shown that in the case of Thomson and 

' Rayleigh scattering, and when scatterers and stars are distributed with ro- 

' tational symmetry, the total polarized flux depends on the inclination, i, of 

the galactic axis to the line of sight according to a simple sin'^ i law. This 
generalises the well known result for pointlike and spherical light sources. By 
using a method based on spherical harmonics, we generalise this law for more 
general mechanisms, and show that to good approximation, the sin^ i law still 
holds for the class of models considered. 

Key words: galaxies - polarization 



1 INTRODUCTION 

The central discs of spiral galaxies are known to contain substantial quantities of dust, 
molecules and free electrons. The scattering of light from stars, which are largely confined to 
the disc and the bulge, by these dust particles, molecules and electrons can broadly explain 
the observed optical intensity and polarization of spiral galaxies. 

Maps of optical intensity of a large number of spiral galaxies have been obtained obser- 
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2 J.F.L. Simmons and E. Audit 

vationally, and polarization maps for a smaller number. Several galaxies have been more or 
less successfully modelled for their intensity and polarization by Monte Carlo simulations 
(( [Wood 1997| )), a technique suited to the optically thick regime. Such studies have yielded 
more information about the dust content of spiral galaxies (( [Byun et al.l99^ )), and in some 
cases have indicated the presence and strength of magnetic fields (( praper et al.l995| ) , 
( ^carrott et al. 1996|) ). Whether the galaxies are optically thin, or optically thick has a cru- 
cial bearing on our understanding of galactic evolution (( Palzetti fc Heckmanl999| )) and star 
formation. Similarly, the detection of magnetic fields would have considerable importance 
for our understanding of galactic evolution. Clearly absorption by dust particles and gas will 
produce an attenuation of light, and hence a decrease in the apparent brightness of galaxies. 
For this reason such calculations have a practical bearing on distance estimation of galaxies. 

Monte Carlo simulations undoubtedly provide a powerful technique for understanding 
the physical processes taking place in the galaxy, but often, because of their dependence 
of a specific choice of geometries and parameter values, can obscure certain fundamental 
properties and relationships that might be obtained through a simple analytic approach. 

In this paper we assume that the distribution of scatterers is optically thin. This as- 
sumption allows one to obtain a number of interesting analytic results for cases where the 
distribution of scatterers and of stars is symmetric, and yields fairly simple expressions for 
the unpolarized and polarized intensity and flux for the more general case in terms of simple 
integrals. In the case of many spiral galaxies the assumption of optical thinness is probably 
not unreasonable (( pCilouris et al.l999|) ) , ( ([Byun et al.l994|) ) but even in those cases where 
it is not expected to hold, the optically thin results can often give a qualitative picture of 
what is happening. It is also possible to give a semi-analytic treatment of the the case where 
the galaxy is optically thick in absorption, but optically thin in scattering. This we shall 
deal with in a future paper. If the spiral galaxy is considered to have a rotational axis of 
symmetry we expect the direction of polarization to be along (or possible perpendicular to) 
direction of axis of symmetry projected perpendicularly to the line of sight. (Of course this 
axis of symmetry would be broken by the presence of spiral arms, but even so would be 
approximately valid.) It has been emphasised (( [Audit fc Simmonsl999 )) that this orienta- 
tion of the total integrated light polarization of galaxies could play and important role in 
studies of the distribution of dark matter from weak lensing, which has the effect of changing 
the orientation of the semi-major axis of the elliptical isophotes of the galaxy, but leaves 
the direction of polarization unchanged. The difference between the direction of the image 
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The optical polarization of spiral galaxies 3 

semi-major axis, and the polarization direction thus gives and indication of the strength of 
the lensing. This potentially would considerably reduce the uncertainty in inferring the mass 
distribution with in the lens compared with the usual weak lensing studies, which take the 
orientation of the source galaxy as unknown. We show in this paper that in this symmet- 
ric and optically thin case for the case of Thomson and Rayleigh scattering obeys a sin^ i 
law, where i is the inclination of the axis of symmetry to the line of sight. This generalizes 



the well known result for Thomson scattering for point light sources sources (Prown et al 



19771) , ( ^immonsl982| ), and for spherical extended sources ( Passinelli et al. 1986|) derived 



in the stellar context. We further show that even for other more realistic scattering phase 
functions this law approximately holds for typical galaxy models, generalising the results of 
( ^immonsl983| ). 

The structure of paper is as follows. In section 2 we introduce a widely accepted paramet- 
ric model for spiral galaxies. In section 3 we give the basic definitions of Stokes intensities 
and fluxes and set down the equations of radiative transfer for the Stokes intensities, and 
specialise these to the case where the source of light is provided by a distribution of stars. 
In section ^ we discuss the optically thin case, and derive expressions for the Stokes inten- 
sities and fluxes for Rayleigh and more general scattering mechanisms, and compare the 
polarization maps with those obtained by Monte Carlo techniques. We go on to outline a 
general and very powerful method for treating the optically thin case for general scattering 
mechanisms that makes use of spherical harmonics and their properties under the rotation 
group. Finally we present our conclusions in section ^. 

2 MODEL GALAXIES 

In this section we outline the model which we have adopted for spiral galaxies. The galaxies 
are taken to be composed of a disc and bulge, which contain stars and scattering dust and 
particles. We take the distribution of stars ( |Jatfe 19831) in the bulge to be given by 



where r is the distance to the galactic centre and is a characteristic radius. Following 
other studies (( |Bianchi et al. 1996| ), ( |Wood fc Jones 1997[ )) we use a value of = 1.0 kpc. 
We have also, for numerical reasons, applied a cut-off of the bulge distribution at a distance 
of 2.5 kpc. 
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4 J.F.L. Simmons and E. Audit 

The distribution of the stars in the disk is given by: 

Pd{R, Z) = pyxp{-R/Rd) exvi-ZjZi), (2) 

where R is the distance to the axis of symmetry of the galaxy and Z the distance to the 
galactic plane. The typical scale length for the disc is taken to be = 4 /cpc with a cut-off 
at 20 k-pc. The typical thickness is Z^ = 0.5 kpc with a cut-off at 2.5 kpc. 

The distribution of the scattering particles (dust or electrons) is also taken to be of the 
same form as equation (|]), i.e. 

n{R, Z) = noexp{—R/ Rg) exp{—Z/ Zg) (3) 

The dust disc is assumed to have the same radial extend as the stellar disc (i.e. Rg = 
R^ = 4 kpc, with a 20 kpc cut-off), but to be thinner: Zg = 0.25 kpc with a cut-off at 
1.25 kpc. The total content of dust is normalized by stipulating that the total optical depth 
of the galaxy along its axis of symmetry. With the form given by equation ^ this is given by 
r = ZgUQcr, where a is the total scattering cross section. In our numerical models we shall 
take T = 0.05. Since we are in the optically thin approximation, results for other optical 
depths can be obtained with a linear scaling. 

3 EQUATION OF RADIATIVE TRANSFER FOR THE STOKES 
PARAMETERS 

In this section we set down the equations of radiative transfer, and apply them to our model 
spiral galaxy. We are essentially interested in the (asymptotic) Stokes intensities and the 
polarized and unpolarized fluxes as measured by a terrestrial observer. 

The brightness, and the degree of linear and circular polarization of a radiation field is 
be described by the Stokes intensities, denoted by /, Q, U, V. I is the unpolarized intensity, 
Q the difference in intensity measured by a polaroid aligned in two perpendicular directions, 
say X and y, U the difference in intensity measured in two perpendicular direction at 45*^ to 
the X and y axis, and V the circular polarization. I,Q,U,V are each functions of position, 
r, and direction n. For convenience we sometimes use the notation Ji = I, I2 = Q, I3 = U, 
and J4 = V, and introduce the 'vector', I = {Ii, I2, 13, h)^ for the Stokes intensities. The 
vectorial Stokes fluxes are defined as the integrated intensities, viz. 

Fi= [ I, n'dn^,. (4) 
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The optical polarization of spiral galaxies 5 

The flux, corresponding to polarization state i, across a surface oriented in direction n will 
then be n ■ Fj. We shall be interested in the flux across a surface in the direction of the 
observer. This is a scalar quantity. We shall sometimes use the obvious notation Fi = Fj, 
F2 = Fq, F3 = and F4 = Fy, and F = {Fi, F2, F3, F^)'^ . To simplify notation, we 
shall implicitly assume that the Stokes intensities can depend on wavelength rather than 



use a lambda subscript. The degree of linear polarization is given by y Fg + F^/Fj, and the 
position angle of the polarization by 1/2 arctanFg/Ffy. 

The direct source of light is supplied by the stars, and this light is scattered by electrons, 
or absorbed, scattered or emitted by molecules and dust. Emission from galactic clouds 
could easily be incorporated into the analysis, as can absorption, but for simplicity we ignore 
these. We assume that the starlight is unpolarized, although this too could be included in 
the present formalism. Throughout this paper we take the galaxy to be optically thin. This 
assumption appears to the valid for most spiral galaxies ( pCilouris et al.l999D ,( Posma et af 
1992|) , although in some cases might break down ( |Scarrott et al. 1996| ). 



Optically thick cases have previously been treated using Monte Carlo techniques ( [Bianchi] 
et al. 1996] ), ( [Wood fc Jones 1997| ), but the purpose of this paper is to treat the simplified, 
though realistic, problem by the simplest techniques. The case where the scatterers are 
electrons or Rayleigh scatterers is largely amenable to simple analytic treatment. In all our 
calculations we assume, fairly realistically, that the stars are unpolarized sources. Usually 
we take a model in which both stars and dust are continuously distributed with rotational 
symmetry about an axis of symmetry of the galaxy. 

The equation of radiative transfer in its full generality may be written 

^ 9I(^n^ + n. V I(r, n, t) = C(r, n, t) - D(r, n, t). (5) 

Equation (H) is of course a set of four partial differential equations for the four Stokes 
intensities I,Q,U, V. We consider only time independent solutions, and thus all quantities 
will be taken to be independent of t. C(r, n) is the energy of the corresponding polarization 
state scattered or emitted into direction n per steradian per unit time per unit volume, and 
D(r, n) is the corresponding energy removed per steradian per unit time per unit volume at 
position r. Generally D(r, n) will be linear in I(r, n). It is convenient to write 



C(r, n) = Cstars(r, n) + Ccmiss(r, n) + Cscatt(r, n) (6) 
where Cgtarslr^n) is the contribution from stellar light sources, which we would expect to 
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6 J.F.L. Simmons and E. Audit 

be isotropic (i.e. independent of n). Ccmiss(i', n) is the contribution from emission processes. 
Stimulated emission would depend linearly on I(r, n) (such a term is easily incorporated 
into the following analysis by simple incorporating it in the term D(r, n)). Thermal emission 
should be isotropic and thus can be incorporated into the stellar term. Here however we shall 
for simplicity ignore Cemiss(i"5 n). Cscatt(i'5n) is the contribution from scattering. 

If we consider a mean stellar luminosity of L, and a stellar number density denoted by 
p(r), then, for unpolarized sources Cstars(r5 n) = (p(r)L/47r, 0, 0, 0). Photons will be scattered 
from all directions into the direction n, and so Cscatt(r, n) will depend on the value of the 
intensity in every direction, n', at r. Thus we are dealing with an integro-differential equation, 
which can, except for exceptional cases, only be solved numerically. It is natural to consider 
the solution of equation (^) along rays (characteristics). The parametric equation for the 
ray passing through some arbitrary point with radius vector fq is r = fq + sn, where s is 
the distance of r from tq along the ray. Introducing the notation r' = fq + s'n, equation (^ 
now takes the form 

^^^ = C(r',n)-D(r',n), (7) 

which has the formal solution 

I(r, n) = I(ro - tn, n) + ^ C(r', n) - D(r', n)rfs', (8) 

where —t is the parameter value at the initial point of integration along the characteristic. 
Allowing t oo, and using the boundary condition I(ro — tn, n) as t oo, equation 
(|) becomes 

l{r,n)= r C,tars(r',n)rfs'+ r Cscatt(r', n) - D(r', n) rfs'. (9) 

J~oo J —oo 

The term D, the extinction, is simply given by D(r, n) = 'n,(r)(Tl(r, n), where a is the 
scattering cross section and n(r) the number density of scatterers. Cscatt is more complicated, 
since it involves an integral of the Stokes intensities over all incoming directions, n'. Moreover, 
the scattering plane for photons scattered into the line of sight n varies with the incident 
direction, and so the total contribution has to be obtained from an appropriate rotation to 
some reference plane. Thus consider an incoming photon from direction n' that is scattered 
by an electron or dust particle at position r into the direction n (see figure ^ . n and n' 
define the scattering plane, the normal to which is given by the vector n x n'. It is natural 
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r 




Figure 1. The incident photon direction n', and scattered photon direction, n, define athe scattering plane. 1' is taken to be 
the normal to the scattering plane, and m' completes the right hand frame {l',m',n'} for the incident photon. The scattered 
photon frame is defined by the right hand frame {1, m, n}, where 1 is taken to be equal to 1'. 

to define a right handed orthonormal basis {!', m', n'} associated with the incoming photon, 
where 1' = n x n'/|n x n'| is the unit normal to the scattering plane and m' = n' x I'. The 
right hand basis {I, m, n}, where 1 = 1' and m = n x 1 is associated with the scattered 
photon. 

We denote by I' the Stokes parameters for incoming photons measured in the frame 
{1', m', n'} . Incoming photons in solid angle dVt-^i scattered into direction n give a contri- 
bution, (il(r,n), per unit volume to the scattered intensities in frame {1, m, n} given by 
dl = o"r2(r)S(n, n')r(r, n')(inn'- This defines the scattering phase function S(n, n'). To in- 
tegrate the different contributions to the scattered intensities from different incoming rays, 
we need to express all the scattered Stokes intensities in the same fixed reference frame. We 
take this fixed frame to be the observer's frame {e^;, e^, e^}, which is chosen such that is 
the direction from the galaxy to the observer, e^. lies in the plane of the axis of symmetry 
of the galaxy and the line of sight, and ey completes the right hand basis. Since in this 
context we are only interested in photons scattered towards the observer, we may without 
loss of generality put n = e^. Thus in order to sum the contributions from different beams 
we need to rotate from basis {1, m, n} to {e^, e^^, e^}, that is through and angle given by 
cos (p = By ■ m. Under such rotations the Stokes parameters transform as 
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(10) 



where 



R 



10 

COS 20 — sin 20 

sin 20 cos 20 

1 



Thus equation(g) becomes 

I(r, n) = Io(r, n) + cr 



ds'n{Y') / R(0)S(n,n')I(r',n')rffin' 



47r 



n{v')a\{v\ n)ds', 



(11) 



(12) 



where Io(r, n) = /fg^ Cstars(r', n)(is' is the intensity at r in the direction n (i.e. when no 
scattering or absorption is present). (Any cyhndrical distribution, p{R,z), of stars, which 
we assume also radiate unpolarized radiation, must necessarily give rise to a source field, 
Io(r, n), that is cylindrically symmetric and unpolarized.) We may easily express the source 
intensities in terms of p{R, z). Indeed, from figure H it can be seen that 



poo 

/ p(r - sn)rfs, 0,0,0). 
Jo 



An Jo 

It is convenient to introduce the normalised stellar density, 

P 



P 



Jpd'^r' 



Equation 13 can then be written 



Io(r,n) = -^(/ p(r-sn)rfs, 0,0,0), 
An Jo 



where Lq is the galaxy's luminosity. Let us define the normalised source intensity 

An 

E(r,n) = — Io(r,n). 



(13) 



(14) 



(15) 



(16) 



S(r, n) is essentially a stellar surface density. 

We shall consider several scattering mechanisms. For Thomson scattering by electrons 
the scattering phase function S takes the form 



^1 + cos^X "'"^ 



16n 



sm X 

2., 1 , 



(17) 





sin"^ X 1 + cos^ X 
2cosx 

2cosxy 
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X 



Figure 2. The intensity in direction n' at rr due to stellar distribution is obtained by integrating the stellar density along the 
line in the direction of n'. The axis of symmetry is Z. The polar and azimuthal angles of n' in the the cylindrical frame are 
given by a and /3 respectively. 

In the case of Rayleigh scattering (molecules and small dust particles) the form of the 
scattering phase function is the same. However the total scattering cross section behaves as 
a oc A""^. 

In many numerical simulations the Henyey-Greenstein form for scattering is used as an 
approximation to a typical mix of galactic dust particles. 
In this case the scattering matrix takes the form 





(p. 


P2 





\ 


1 


P2 


Pi 


















-Pa 











P.) 



(18) 



where 



and P2 is given by 

J-, Pmax -1 . 9 

Pi 1 + 

P3 and P4 are irrelevant for our study, since we are only concerned with unpolarized 
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incident light. The single parameter which represents the mean value of the cosine of the 
scattering angle, determines how peaked in the forward direction the scattering is. Pmax is 
the peak polarization (i.e. the polarization at 90°). For the scattering of optical light on 
dust, g and Pmax are both of the order of 0.5 ( [White 1979| ). We use these values in all our 
numerical calculations. Throughout, we have assumed that there is no absorption and that 
attenuation is purely due to scattering, corresponding to the case of albedo 1. 



4 THE OPTICALLY THIN APPROXIMATION 

If the galaxy is optically thin, in the sense that / nadl through the galaxy is less than one 
in all directions, then the first term on the right hand side of equation (|1^) will dominate. 
The optically thin approximation is obtained from first order iteration, in which I(r, n') in 
the integral over solid angles in the second term, and 10(1", n') in the third term are replaced 
by and Io(r, n') and Io(r, n) respectively. Let us introduce the notation dTz = an{r)dz. Thus 
in the optically thin approximation the asymptotic Stokes intensities in the direction of the 
observer, I{x, y, 00, e^), are given by 

rr{x,y) r ^ ^ 

I{x,y,oo,ez) =lo{x,y,oo,ez) + dr^ R{4))S{ez,n)lo{r,n)dQn' 

Jo J 4. ^^g) 

- / Io(r,eJ (ir^. 
Jo 

where T{x,y) is the optical depth through the galaxy at the field point {x,y). 

Writing equation (^) in component form, and introducing the elements Sij of the scat- 
tering phase function, S, we obtain 



fr(x,y) r ^ ^ 

I{x, y, 00, e^) = Io{x,y, 00) + dr^ sii{n,n)Io{r,n)dQn' 

Jo J4n 

fr{x,y) 

- / Io{r,ez)dTz 
Jo 

Q{x,y,oo,ez) = J dr^ rffin'/o(r, n')(s2i(n, n) cos2<^ - S3i(n, n) sin2<^) (21) 
U{x,y,oo,ez) = J ^Tz J (ifin'/o(r, n')(s2i(n, n') sin20 + 531(11, n') cos20) (22) 
V{x,y,oo,ez) = / dr^ / c/iln//o(r, n')s4i(n, n) (23) 

J J4n 



(20) 



and 
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where Jq is given in terms of the stellar number density by equation (|I3|). Introducing the 
notation 



W = Q + iU (24) 

we obtain 

W{x,y,oo,e-^) = dn (if2n'/o(r, n') expi20((s2i(n, n') + is3i(n, n')). (25) 

J JAtt 

In the case of spherical scattering, where the angular dependence of Sij is only in the 
scattering angle, 6', one can expand Sn {s^i = in this case) in terms of Legendre polynomials 
P;, and S21 and S31 in terms of associated Legendre polynomials, Pim with m = 2. This leads 
to simple expressions for the asymptotic Stokes intensities in terms of the line integrals of 
the multipole moments of the Io(r, n), with m = for the I and V and m = 2 for Z. This is 
similar to the results obtained by ( ^immonsl982|) and ( ^immonsl983| ) for the point source 



case. For Thomson and Rayleigh scattering further simplification results from the fact that 
only terms with / = 0, 2 occur in the expansion of the phase functions Sij. This has a crucial 
importance when we come to calculate the flux. We shall not pursue this line of reasoning 
here, but rather adopt the simpler approach adapted to the case of Thomson and Rayleigh 
scattering. However, we outline the method using spherical harmonics section |4.3| . 

4.1 Polarization map 



Using equations ^0H23|, it is possible to compute polarization maps of the galaxy. To carry 
out this we have divided the galaxy field into 100 x 100 pixels, and for each each pixel 
we have integrated along the line of sight. The results our shown in figure ( [4.1|) , and can 
be seen to be comparable with those obtained by Wood ([Wood 1997| ), who used Monte 
Carlo techniques. Of course the integration is in our case extremely quick. The degree of 
polarization is much higher for a galaxy viewed edge on than for the same galaxy viewed 
face on, and this is essentially due to the respective optical depths. One should note that 
for Henyey-Greenstein scattering the upper and lower half of the galaxy polarization is no 
longer symmetric. This is due to the functional form of the scattering, which is now peaked 
in the forward direction. 

One can obtain the polarized flux from the polarization intensity maps by integrating 
the contributions to the polarization over the field of view. However, as we show in the 



sections [4.2| and [4.3| , analytic expressions for the polarized flux in the case of Thomson 
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Figure 3. Polarisation maps for Thomson scattering (bottom plots) and Henyey— Greenstein scattering (top plots). The in- 
clinations of the galaxy are 87°, 60° and 18° for the left, central and right plots respectively. The degree of polarization is 
proportional to the length of the lines, the maximum polarization being 0.8% , 0.5% and 0.2% for the top row and 1.8%, 1.1% 
and 0.4% for the bottom row. The axis a labelled in kpc. 

and Rayleigh scattering can be found, and approximate expressions for general spherical 
scattering mechanisms. 

4.2 Polarized flux for Thomson and Rayleigh scattering 

The flux at the earth is obtained by integrating the asymptotic intensities, i.e. y, oo, e^^), 
for large s over all solid angles as seen by the observer. Noting that this solid angle is just 
given by dxdy/R\ it is clear that integration of equation (|l^) over solid angles gives an 
expression for the flux in terms of volume integrals and the direct flux, Fo(c>o). 

F = ^ J I{x,y, oo, e^)dxdy = Fq + J dVan{v) J R(0)S(e^, n')Io(r, n')dl]„/ 

~ J (Tn(r)Io(r, e^)rfV. ^^26) 
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X in kpc 



Figure 4. Polarization profile along the major axis of a galaxy for Thomson scattering for an inclination 87°. For Henyey— 
Greenstein the profile is very similar but the overall degree of polarization lower. The cut-off at ±20fcpc corresponds to the 
cut-off we have applied in the disc distribution, but the intensity at this distance from the galactic centre is already extremely 
low 



From equation ([TsD it follows that 



Fo(oo) = ^(1,0,0,0) (27) 

E 



When calculating the degree of polarization to first order in optical depth, we need 
only normalize Fu and Fq by the direct unpolarized flux, ignoring the contribution to the 
unpolarized flux from scattering. 

Let us now assume that the galaxy has axial symmetry. Thus the distribution of stars 
and the density distribution of scatterers are taken to be cylindrically symmetric about this 
axis. The axis of symmetry is inclined at an angle of i to the line of sight. It is useful to 
introduce cylindrical coordinates with origin at the centre of the galaxy, {i?, $, Z} , with 
the Z direction oriented along the axis of symmetry of the galaxy. (See figure ^]2| .) 

It is also convenient to write n = no n, where no is the central density of scatterers, and 
h thus a dimensionless density. We introduce further a length scale Zg and the dimensionless 
variable Z' = Z/Zg, which will enable us to express the Stokes fluxes in terms of the optical 
depth. To = n^aZg. This is indeed the vertical optical depth for the models of section H. We 



shall also work with the normalised intensity, S(r, n), given by equation [16 



Because of the symmetry, the density of scatterers, n, is only be a function of R and 
Z. Indeed this is the form assumed in the models of section |[ Similarly the unpolarized 
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Figure 5. The axis of symmetry, given by Z, is inclined at angle i to the line of sight. The observer's frame is {ei^,ey,ez}. 
OX, OZ and are taken to lie in the same plane. 

source Stokes intensity, S(r, n), is a function of R and Z, and the direction cosines, a, (3 of 
n', expressed in terms of the associated coordinated basis vectors {e/j, e$, e^}. 

We shall now prove that under these assumptions, and with the additional assumption 
that the dust is optically thin, Fq depends only on sin^ and Fu = 0, generalizing the result 
obtained ([Brown et al. 197^) for point sources. The behaviour of the degree of polarization 



will differ slightly from this, since the normalization factor, Fj, also depends weakly on the 
inclination. However, to first order in the optical depth the result will still hold. We also 
obtain an explicit form for the polarized flux in terms of density moments. Equation (^6|) 
written in component form now becomes 



z,a,f]) sin a da dp 
-J n{R,Z)aIo{R, Z,a,P)d\^28) 

Fq = — %r [ h[R, Z)RdRdZ'd(j) [ sin^x S(i?, Z, a, /5) cos20 smadadp 
4:71 lovrrip J J 



(29) 



and 



Lg 3ro 
Uj = TT / n 



{R, Z)RdRdZ' d(j) j sin^ x ^o{R, a, j3) sin 20 sin a da dj3. 

(30) 
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To evaluate F/, Fq and Fu we shall need to work out the scattering angle x (between n = 
and n') and the rotation angle between the scattering plane and the observer's y axis in 
terms of $,a;,/3. (As we shall see, although we do not require an explicit expression for the 
scattering angle in terms of for the calculation of Fq and Fjj, for the evaluation of 

Fi we shall require it.) 

In order to do this we need to express the basis {e^j, e$, e^} in terms of {e^., e^,, e^}. Let 
us introduce the convenient shorthand notation for cosine, cos A = c^, and sine, sin A = sa- 



Cq> Sq, 
— S$ Cq, 

V V 



0\ ( C,; (e.y 



1 

-Si Ci 



■y 



(31) 



The first matrix on the right hand side represents a transformation from {eij,e$,e2} to 
{ex, ey, e^}, and the second from {ex, ey, e^} to {e^^, e^, e^}. 
Equation (^Tj) simplifies to 



e$ 
\ezj 



~ S ^ C^i S ^S'l 

Si Ci 



■y 



(32) 



Writing 



n 



SaCp SaSp Cc 



(33) 



and using equation (0) we obtain 

II ~ (-^aC^CciiCj S(y^SpS^Ci -\- C(y^Si)Gx 

+ {SaCpSq. + SaSpCq,)ey (34) 
+ (SaC/3C$Si + + CaCi)ez 

Noting that coscj) = (n' x e^) ■ Gy/ln' x e^l and sin0 = (n' x e^) ■ e^/|n' x e^l we immediately 
obtain from equation ( [5^ ) 

Ca^i ~l~ SaCpC^Ci S^CiSaS g 



COS( 



and 



sm I 



(35) 



(36) 



The scattering angle, x, is given by cosx = (n' ■ e^), and again from equation ( [3^ ) we obtain 

— CaCi S^yCpCqSi ~1~ Sf^SqyS ffSi. (^'^) 
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16 J.F.L. Simmons and E. Audit 

Noting that cos 20 = cos^ (p — sin^ 0, and sin 20 = 2 cos sin and using eqs. (|35|) and 



we see that in the integrand of equations ^ and ^ the terms sin^ x cos 20 and sin^ x sin 20 
simphfy by virtue of the cancellation of the sin^ x terms, and yield respectively 

and 



"^X"^20 



-2(CaSi + Sa^CpC^Ci - S$QSaS/3)(SaC/3S$ + SaS/jC^). (39) 



Y h{R,Z)RdRdZ' J /o(i?, Z, a, /3)(3c^ - 1) Se^dpda. (40) 



Upon carrying out the integration over $ all terms in sin$, cos$, and sin 2$ disappear, 

since both n and / are independent of $. A calculation shows that only terms with the 

factor sin^ i remain in the final expression for Fq, and Fu = 0. The latter implies that the 

polarization is in the y direction or in the x direction, i.e. along or perpendicular to plane 

defined by the axis of symmetry and the line of sight. Explicitly we obtain 
_ Lg t sin^ i 
~ 47r 16i?| 

Equation ( ^OD has a simple interpretation: the polarized flux is proportional to the dipole 
moment of the radiation field weighted over the density distribution of the scatterers. 

A similar calculation can be carried out for the unpolarized flux. Substituting in equation 
(p8|) the expression for cosx given in equation (0), we obtain 

Fi = Fo(oo) + ^^^^ / n{R, Z)R dR dZ' J /o(i?, Z, a, (3){2 + Icl^^ - sls^c^ii) Sad(3da 

- J n{R,Z)aIoiR, Z,a,P)0^ 
The degree of linear polarization, p, is given, to first order in the optical depth, as 
FqI Fq{oo). From equations (^ODand (|27|) we obtain 



Stttq sin i 



J RdRdZ'n{R, Z) J Z, a, p){3cl - 1) s^dpda (42) 



where Z, a, (3) is given by equation |T6. 



4.3 Approximate expressions for polarized flux using spherical harmonics 

For more general mechanisms than Thomson and Rayleigh scattering the analytic methods 



used in section are not very helpful. The approach used there worked because of the 
cancellation of the factor sin^ x the integrand of the polarized fluxes owing to the form of 
the scattering phase function. For most phase functions this will not be the case. 
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The optical polarization of spiral galaxies 17 

Spherical harmonics, Yim{0,(j)), provide us with a very powerful method of dealing with 
the general case. On the one hand they provide a complete set of orthogonal functions on the 
sphere, and on the other for each value of / the 2/ + 1 function Yim{0, <f)) provide an irreducible 
representation of the rotation group. We shall largely follow the notation of Messiah (1961), 
and take the spherical harmonics to be defined by 

YimiOA) = c{l,m) Pr{cose)e'"^^ where c(/,m) = [i^^±M_^]i/2 (43) 

and are the associated Legendre polynomials. c{l, m) is a normalisation constant to 
ensure Yim{0,4') are orthonormal. 

Thus suppose that we have two coordinate bases, {e^;, 6^,62} and {e^,e^,e'^}. Suppose 
further that the rotation taking {e^.. By, e^} into {e^, e^, e'^} is described by the Euler angles 
A, B, C. We then have the relation 

m'=l 

Yira{e\ct>')= E Yu^,{e,cl>)R^!^,^{A,B,C). (44) 

m'=-l 

The matrices R*-'-* form an irreducible representation of the rotation group. The elements of 
R^'^ take the form 

C) = e-"^-^r«„,(5)e— (45) 
where r'f^,^{B) can be calculated from the Wigner formula (see Messiah 1961), 



r. 



mm 



(I) ^ . ./(/ + m)\{l - m)\{l + m')\{l - m'). , ^ , 

('j V ^ c^l+m~m —2t^2t—m+m 



y^/ V ^ ' ^ / \ • ' \ ' c^l+m~m'-2t 

Y {I + m - ty.{l - m' - ty.t\{t - m + m'y ' 

(46) 

where ^ = cos-B/2 and 77 = sin 5/2. 

It is convenient for our purposes to work entirely in the observer's reference frame (see 
figure ( [4.2|) . Our two frames are thus {e^;, e^,, e^}, the observer's frame, and {e'^,e'y,e'^} = 
{eij,e$,e^}, the frame attached to the galaxy's cylindrical coordinates. With this notation 
we replace 9' by a, and 0' by /3. To rotate the coordinate basis {e^;, By, e^} into {e/j, e$, e^} 
we have to 

(i) rotate about y-axis though i and then 

(ii) rotate about Z-axis though $. Thus A = 0, B = i, and C = $. 

Let us assume that the source Stokes intensity has cylindrical symmetry, i.e. /q = 
Iq{R, Z,a, P). Let us work in terms of the normalised surface density E(i?, Z, a,/3) = 
4:7iIq{R, Z, a, 13) /Lq as defined in equation |16|. Expanding E(i?, Z, a, /3) in terms of spherical 
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harmonics 

l=oo m=l 

E{R, Z,a,P)=J2 E ^imiR, Z)Yim{c^, P) (47) 

1=0 m=-l 

where of course aim do not depend on Indeed because of the orthogonahty of the spherical 
harmonics we may write 

aiUR,Z)=f S(i?,Z,a,/3)r;^(a,/3)da (48) 

Jin 

Now expressing Yim{a, P) in terms ofYim{0, 0) via equation^ and substituting into equation 
i7| , we obtain 

/=oo m,=l m'=l 

E{R,Z,a,P) = J2 ^irn{R,Z)Yim'{9,m^J;!'miO,t,^). (49) 

/=0 m=—l m'=—l 

Now from equation representation, i?£'',^(0, i, $) = T£/^(i) exp —im^. 

The angle ( denoted by (p in figure ^) between the scattering plane of the incoming photon 
in direction n, and the observer's x — z plane, is = — tt/2, and the scattering angle, x 
in the main text, is x = ^• 



The complex intensity, W, is given by equation 25 . Note that for spherical scattering 
mechanisms we need only consider S21 and moreover S21 = S2i{6). To obtain the polarized 
fiux, expression (p5|) for the complex intensity, W = Q + iU, has to be integrated over 
all solid angles at the observer. As in section ( |4.2| ) we multiply by the element of solid 
angle, dxdy/R\, thus transforming the expression (^) into a volume integral, where we 
have written dTz = auQudz . 

As before, introducing the dimensionless variable Z' = Z/Zg, and the central optical 
depth, dro, we thus obtain the expression 

FQ + iFu = / h{R,Z)RdRdZ'd<^> [ S2i{9)E{R,Z,a,p)e'^''' sin 9d9d(l). 

(50) 

Substituting equation (|49|) into equation (|^) we obtain 

1=00 m=l m'=l 



-TO + - 



AnR 



/n t III, (, III, 1, 

h{R,Z)RdR dZ'd^ / S2i{9) E E E 

i=0 m=—l m'=—l 



(51) 

aiUR,min^'{0Ay3j^)e-'"'''e'''^sm9d9d(f) 
Integration over $ eliminates all terms for which m 7^ 0, and introduces a factor 2tt for 
m = 0. With the definition 

Aim = J n{R, Z) aim{R, Z) R dR dZ' (52) 
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equation (|5T| ) reduces to 

T ^ l=oo m'=l „ 

Fq + iFu = E E / S2i{e)Yim'ie,^)r^J^,,{^e'"^smeded<p. (53) 

Clearly, only the term m' = —2 can contribute, and integral ( |5^ ) reduces to 

T l = 00 

Fq + = E ^iAm2 (21) (54) 

^-n-E 1=2 

where 

5z2 (21) = 27r c(/,-2) 1 S2i{e)Pi2{cose)siYiede. (55) 

Since all the terms on the right hand side equation (|5^ are real, it follows that Fjj is zero, 
meaning that the polarization is either directed along the axis of symmetry, or perpendicular 
to it. For Thomson scattering the only non-zero term in the multipole expansion of the phase 
function is 5*22 (21)) and we obtain the result proved in ^]2|, since ri^2o(^) is proportional to 
sin^ i. On the other hand, if the scattering phase function has higher multipoles, then the 
dependence of the polarization on the inclination angle will have higher terms in sini and 
cosi. 

To first order in optical depth, the degree of polarization is given by, p = |Fq|/Fo(oo). 
Noting that -Fo(oo) = Lg/^t^R'e "we obtain for the degree of polarization 

(=oo 

p = 2TxroY.Aior%S)Si2 (56) 
1=2 

A similar calculation can be carried out for F/. 

4.4 Dependence of polarized flux on inclination 

From the polarization maps presented in section (f4.1|) , it is possible to calculate the inte- 
grated degree of polarization over the entire galaxy. This must be a function of the inclina- 
tion. The corresponding plots of polarization (calculated in this way) against inclination are 
given in figure ([4.4|) . For Thomson scattering one obtains the sin^ i law, as expected from the 
analytic results of section |4.2| . On the other hand, for Henyey-Greenstein scattering, there 
is a noticeable deviation from this law. 

The dependence of the total polarization on inclination can also be calculated using the 
expansion in spherical harmonics outlined in the previous section. The degree of polarization 
is given by equation |5^. 

AiQ is an integral over the distribution of stars, dust and gas, and does not depend 
on the scattering mechanisms or inclination. r^%Q is a function of inclination that is easy to 
calculate. 5*^2 (21) is determined from the scattering function from the simple integral given by 
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(0 

^-20 




•S';2 (21) (thorn.) 


Si2 (21) (HG) 


1=2 




2.21 10^2 


^Stt/IO 


0.373 


1=3 


^sm2(i)cos(i) 


-4.16 10-3 


0. 


0.174 


1=4 


^sin2(i)(5 + 7cos(2j)) 


-1.72 10^2 


0. 


0.043 



Table 1. 




Figure 6. Degree of polarization as a function of i. The empty squares show the numerical values for Thomson scattering and 
the full squares for Henyey— Greenstein scattering. The si'n? law is, as was expected, recovered for the Thomson scattering (solid 
line) . For the Henyey-Greenstein scattering there is a noticeable deviation from this law. The short-dashed line corresponds 
to the 1=2 approximation, and the long-dashed line to the 1=4 approximation for Henyey-Greenstein scattering. 



equation |55|.Results for Z = 2, 3 and 4 are shown in the table |l]. For Thomson scattering only 
/ = 2 is non-zero. For Henyey-Greenstein higher order multipoles are not zero, nevertheless, 
as is illustrated in figure ( ^.41 ), even with inclusion only of the / = 2 term a very good fit 
is found. Including terms up to / = 4, one obtains an excellent fit to the numerical results, 
with a disagreement of only a few percent. 

Evidently this method can test an arbitrarily large number of different scattering mech- 
anisms and geometries extremely quickly, since the total polarization is given by a simple 
sum of products of easily calculated terms. Thus to see the effect of changing the scattering 
mechanism, it is sufficient to simply recalculate the corresponding coefficients Sn (21) • 
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5 DISCUSSION AND CONCLUSIONS 

The assumption that the spiral galaxy is optical thin allows us to easily calculate numerically 
the spatially resolved degree of polarization and unpolarized intensity of starlight scattered 
by dust, electrons and gas. Our numerical results agree qualitatively with the observations for 
spiral galaxies, and for vertical optical depths of around 0.05 give maximum polarization of 
around 1% depending on the dominant scattering mechanism and inclination of the galaxy. 
For the same optical depth, Thomson and Rayleigh scattering are more efficient in producing 
polarization than the Henyey/Greenstein phase function. The latter produces an asymmetry 
in the polarization about the semi-major axis of the elliptical image of the galaxy. 

Although we should expect the optical thin assumption might break down for high 
inclination galaxies, which are indeed the ones discussed by ( Praper et al.l995|) and ( ^carrott 



3t al. 199^) , and where the optical depth along the galactic plane is greater than 1, our results 



appear to agree well with the Monte Carlo calculations of a number of authors. It would be 
interesting to make a detailed comparison of the two methods. 

The advantages of the optical thin assumption is undoubtedly its efficiency: the computer 
time needed for calculation of maps considerably less than Monte Carlo treatment. Although 
we have dealt only with Thomson and Henyey-Greenstein type scattering functions, the 
inclusion of mixtures of scatterers etc. would be very straightforward, and involve only 
slightly more computer time. It would is a straightforward extension to treat a case where 
the galaxy is optically thick in absorption, but optically thin in scattering. However in this 
case we should not expect the same inclination law to hold for the polarization. We shall 
deal with this case in a future paper. There have been only few observational studies of the 
polarization of the integrated flux from spiral galaxies. We find that in the optically thin 
regime, for Thomson or Rayleigh scattering, the dependence of the degree of polarization on 
inclination of the galactic axis to the line of sight, has a simple siv? i form if we assume 
axial symmetry in both the distribution of scatterers and stars. Again, for vertical optical 
depths of r ^ 0.05, the total integrated polarization reaches about 1% for higher inclinations. 

Thus the suggestion that optical polarization could be used in the study of weak lensing 
( [Audit fc Simmonsl999| ) is further supported by this study. Although the observation of such 



levels of polarization would be difficult in such high redshift galaxies, the more important 
determination of the direction of polarization might be feasible. 
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Such a technique, in determining the orientation of the source galaxy, would provide 
considerably more precision in determining the mass distribution in weak lensing studies. 

For more complicated scattering mechanisms one can write down an expression for the 
total integrated polarization in terms of a spherical harmonic expansion expansion. For the 
class of galaxy models we have taken, it appears that the convergence of this expansion is 
very rapid, and that to a good approximation the polarization is given by the first term of 
the expansion, / = 2. This indicates that the sin^ i law should be more generally applicable. 
Our formulation easily allows the inclusion of other mechanisms (absorption and dichroism 
etc.). 

Although we have not investigated the the joint distribution of flux and polarized flux for 
a catalogue of spiral galaxies, it would be interesting to do so. In the optically thin regime 
we would expect a correlation between the two, whereas in the optically thick regime the 
flux should be independent of inclination. Thus this could provide a test for optical thickness 
for a class of spiral galaxies. The form of the joint distribution could also indicate whether 
scattering or other mechanisms just as alignment of grains by magnetic flelds is responsible 
for polarization in a class of spiral galaxies. 

Finally we would like to point out the possible use of polarization in a Faber-Jackson type 
distance estimators arising from the possible correlation between polarization and absolute 
magnitude of galaxies, which could be used to reflne the distance scale for galaxies. 
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